Monovalent impurities on graphene: midgap states and migration barriers 
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Monovalent impurities on graphene can be divided into ionically and covalently bond impurities. 
The covalent impurities cause universal midgap states as the carbon atom next to the impurity is 
effectively decoupled from the graphene ir bands. The electronic structure of graphene suppresses 
migration of these impurities and making the universal midgap very stable. This effect is strongest 
for neutral covalently bond impurities. The ionically bond impurities have migration barriers of 
typically less than 0.1 eV. An asymmetry between anions and cations regarding their adsorption 
sites and topology of their potential energy landscape is predicted. 
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The discovery of graphene [l[ and its remarkable elec- 
tronic properties 0,13, m initiated great research interest 
in this material. Particularly prospective for applications 
is the extraordinarily high charge carrier mobility fi in 
graphene P, @, Q . In combination with a very high Fermi 
velocity v ~ 10 6 m/s this makes micron mean free paths 
routinely achievable. 

Away from the neutrality point, the conductivity 
of graphene is weakly temperature dependent and ap- 
proximately proportional to the carrier concentration n 
0, @] ■ The mechanism limiting the electron mobility in 
graphene is still under debate. Charged impurities are 
probably the simplest and thus the most natural candi- 
date [IGSELli- However, room-temperature experiments 
with gaseous adsorbates such as NO2 have showed only a 
weak dependence of // on charged impurity concentration 
Furthermore, recent experiments [l3| did not find 
any significant dependence of n on immersing graphene 
devices in high-K media such as ethanol and water (di- 
electric constants k « 25 and 80, respectively) but this 
disagrees with another report L4| in which a few mono- 
layers of ice increased n in graphene by ~ 30%. Because 
of the experimental controversy, alternative mechanisms 
such as scattering on frozen ripples [loT ] and resonant im- 
purities [3, 17 1 were discussed. However, both these 
mechanisms also have some weak points, rather theo- 
retical than experimental. There is still no real theory 
which would explain why the ripple structure becomes 
quenched and thus almost temperature independent. As 
for the resonant scatterers, in general, the closeness of 
the impurity quasilocal states to the neutrality point nec- 
essary to make this mechanism efficient looks just acci- 
dental. To discuss this as the main scattering mechanism 
one needs, at least, to clarify the mechanism which makes 
the resonant scatterers typical for graphene. Thus, the 
physics of charge carrier mobility in graphene, crucially 
important for most of potential applications, is not clar- 
ified yet. Also, while impurities appear as undesirable 
residua from the graphene production process, chemical 
functionalization of graphene relics on impurities for con- 
trolling its electronic properties as demonstrated recently 



for hydrogenated graphene (graphane) [18| . 

For judging which impurities might determine electron 
scattering in graphene and for optimizing chemical func- 
tionalization, the mechanisms determining the impurity 
mobility need to be known. In this letter, we consider 
monovalent adsorbates and show that these can be di- 
vided into two separate groups regarding the bonding 
mechanism: ionically and covalently bond impurities. To 
this end, we present ab-initio calculations on H, Li, Na, 
K, Cs, F, CI, Br, I, CH 3 and OH adsorbates on graphene. 
For these systems the electronic structure and migra- 
tion barriers are analyzed. The covalently bond impu- 
rities cause a characteristic midgap state derived from 
the graphene electrons. This state turns out to be very 
stable, as graphene's conjugated n bonds enhance the 
migration barriers of neutral covalently bond impurities. 

For a first principles description of the graphene ad- 
sorbate systems we performed density functional calcu- 
lations within the generalized gradient approximation 
on 4 x 4 graphene supercells containing 
one impurity. The Vienna Ab Initio Simulation Package 
(VASP) [2l[ with the projector augmented wave (PAW) 
[22l [23| basis sets has been used for solving the resulting 
Kohn-Sham equations. In this way we obtained relaxed 
structures for the graphene adsorbate systems, total en- 
ergies, and orbitally resolved local density of electronic 
states (LDOS). 

The local electronic structure of graphene in the vicin- 
ity of adsorbates (Fig. QJ can be grouped into two classes. 
The LDOS in the vicinity of adsorbates like Li or CI ex- 
hibits a sharp resonance close to Fermi level which is 
almost entirely localized at the impurity. Besides this 
peak, the LDOS at the nearest neighbor and at the next 
nearest neighbor of the impurity exhibits the pseudogap 
characteristic for graphene. This is qualitatively differ- 
ent for the second group of impurities (Fig. [T] c) and 
d)). H and F adatoms cause a midgap state characteris- 
tic for Dirac fermions: With the bonding partner of the 
impurity in sublattice A the impurity state is localized 
in sublattice B and at the impurity atom. 

Every stable atomic configuration under investigation 
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Figure 1: (Color online) LDOS in different graphene adsor- 
bate systems, (a-b) ionically bond impurities, (c-d) covalently 
bond impurities, a) Graphene + Li, b) Graphene + CI, c) 
Graphene + H d) Graphene + F. For the impurity's bonding 
partner in graphene and its nearest neighbor the p z projected 
LDOS is shown. The valence electron LDOS at the impurity 
site for ionically bond impurities is depicted in (e) and for 
covalently bond impurities in (f). In (a-d) and (f) the Fermi 
level is at E — 0; in (e) the Dirac point is at E — 0. 



can be strictly grouped either into the class of strongly 
or weakly hybridized impurities, as can be seen from Fig. 
[T] e)-f ) : Ionically bond impurities give rise to a sharp 
acceptor (donor) level below (above) the Dirac point at 
Ejy = 0. The LDOS of covalently bond impurities is 
much broader and exhibits characteristic resonances far 
below the Fermi level (between — 10 eV and — 4eV) as 
well as a midgap state at the Fermi level. 

As regards electron scattering this midgap state is 
mainly independent of the particular type of covalent im- 
purity. The supercell band structures for H, CH3, OH, 
and F covalently bond to graphene are shown in Fig. [5] 
The band structure of graphene with adsorbed H and 
CH3 as well as those for graphene with F and OH adsor- 
bates coincide close to the Fermi level despite the differ- 
ent internal structure of the adsorbates. The coupling of 
the midgap state and the graphene bands can be quan- 
tified in an effective impurity model H = Hd + H lmp , 
where the unperturbed graphene bands are described by 
Hb = ^2 k e(k)d' k dk and the perturbation by 

-ffimp = Cimp^C + ^ V k^d k + h.C. 
k 

Here, the index k = (fc, v) denotes crystal momentum 




Figure 2: (Color online) Supercell band structure of a) H and 
CH3 on graphene as well as b) F and OH bond covalently to 
graphene. Bands coincide close to the Fermi level (E = 0) 
despite the different internal structure of the impurities. 



Table I: Impurity energy and hybridization for the effective 
impurity model of the midgap state for different impurities. 
All impurities are placed on top of a C atom, which is at the 
total energy minimum for the covalently bond impurities and 
the anions CI / Br but not for the cations Li and Na. 





£imp (eV) 


\V\ (eV) 




£im P (eV) 


\V\ (eV) 


H 


-0.03 


1.38 


Li 


1.17 


0.22 


CH 3 


-0.11 


1.40 


Na 


0.93 


0.13 


OH 


-0.70 


1.30 


CI 


-0.79 


0.43 


F 


-0.67 


1.30 


Br 


-0.73 


0.19 



k and band number v = ±. e(k) is the unperturbed 
graphene dispersion. The effective impurity is character- 
ized by its energy e imp and its hybridization 14 with the 
graphene bands. In a supercell calculation at the back- 
folded Dirac points k — , this model simplifies to 




(1) 



where the zero block stems from the graphene bands at 
the Dirac point and A (B) are the components of Vk in 
the two different sublattices A (B). This allows to derive 
the coupling strengths | V\ — \/\A\ 2 + \B\ 2 and ei mp from 
the DFT energies of the three bands closest to the Fermi 
level in the supercell calculation as shown in table |U 

The coupling \V\ is mainly independent of the inter- 
nal structure of the covalently bond impurities. Hence, 
this midgap state appears as a universal feature of all 
monovalent impurities which are strongly bond to one of 
graphene's carbon atoms. 
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Table II: Migration barriers AE and minimum energy site for 
ionically and covalently bond impurities. 





site 


AE (eV) 




site 


AE (eV) 




site 


AE (eV) 


H 


t 


1.01 


Li 


h 


0.31 


CI 


t,b 


< 0.005 


CH 3 


t 


0.63 


Na 


h 


0.09 


Br 


t,b 


< 0.005 


OH 


t 


0.53 


K 


h 


0.06 


I 


t,b 


< 0.005 


F 


t 


0.29 


Cs 


h 


0.04 









Bonding H-atoms to graphene and related electron 
scattering has been analyzed in 24, 25, [26|. For the 
bonding partner of H, the 7r bond to its nearest carbon 
neighbors is broken and a a bond with the H ad-atom 
is formed. The carbon bonding partner of H atoms has 
been found to be decoupled from the graphene 7r-electron 
system and the resulting local imbalance between the 
number of atoms belong to each of the two sublattices 
causes a midgap state. The band structures from Fig. [2] 
and the coupling constants from Table |l] show that same 
mechanism is effective for all monovalent covalently bond 
impurities on graphene. 

In the following we show, that the creation of this 
midgap state by an impurity covalently bond to one 
carbon atom enhances migration barriers for covalently 
bond monovalent impurities in graphene. A compari- 
son to ionically bond impurities is given. To find migra- 
tion barriers for ionically bond impurities it is sufficient 
to perform structural relaxations with the impurities in 
three different high symmetry adsorption sites: on top 
of a C-atom (t-site), in the middle of a hexagon (h-site) 
and above the middle of a nearest neighbor C-C bond (b- 
site). The covalent impurities cause strong distortion of 
the nearby bonds and require the minimum energy paths 
to be calculated using the nudged elastic band method 
27|. 

In agreement with [28| . we find the energy minimum 
for the alkali cations at the h-sites and barriers as shown 
in table HIl The barriers decrease with cation size and are 
all (except for the special case of Li) below 0.1 eV. The 
potential energy landscape for the cations consists of dips 
in the center of the hexagons bordered by a hexagonal 
net of banks. Within this net spanned by the nearest 
neighbor carbon bonds, the variation of potential energy 
is by a factor of more than 5 smaller than between the 
h-site and the t/b-sites. 

This landscape is reversed for the anions: Having their 
energy minima on the net and maxima in the center of 
the hexagons, the anions can freely move on the graphene 
sheets. The fact that the height of the impurity above the 
sheet is always minimized in center of the hexagon, would 
result in the energy minimum being in the center of the 
hexagon for all ionically bond impurities if atomic scale 
inhomogeneities in the screening charge of the impurities 
were negligible. The anions preferring the t- and b-sites 



over the h-sites shows that inhomogeneities in the screen- 
ing charge corrugate the potential energy landscape of 
the ions on the order some 10 meV. 

While this can be quite significant in the case of purely 
ionically bond impurities our calculations show that the 
potential energy landscape for covalently bond impurities 
is by an order of magnitude more corrugated. We find mi- 
gration barriers between 0.29 eV for F and 1.01 eV for H. 
Notably, F is the impurity with the highest absolute bind- 
ing energy (1.99 cV) of all impurities considered, here, 
but it has the smallest migration barrier within the group 
of covalent impurities. For F and OH we find the saddle 
point energy of the transition path significantly below 
the desorption energy, which is 0.91 eV for OH. This is in 
strong contrast to H and CH3: For H the energy of the 
saddle point state is only 4 meV below the desorption 
barrier and moving a CH3-group from one carbon atom 
to its nearest neighbor requires even overcoming the des- 
orption barrier of 0.63 eV. No saddle point configuration 
with the CH3-group in the middle of two neighboring C- 
atoms except for the CH 3 being fully desorbed from the 
graphene sheet could be found. 
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Figure 3: (Color online) Electronic structure in the transition 
state, a) p z -LDOS at a carbon atom next to the impurity. 
Lower part: Band structures of graphene supercell with b) H 
and c) F impurities in transition state configuration. Contri- 
butions from the impurity atoms are marked as fatbands. 



For F, OH, and H the nature of chemical binding in 
the saddle point configuration can be understood from 
the supercell electronic properties shown in Fig. [3] The 
LDOS at the carbon atoms next to the F and OH im- 
purities is very similar to the LDOS in the vicinity of 
ionically bond impurities like CI or Br (see Fig. [TJa). 
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This is in contrast to the case of H, where in addition 
to a resonance at 2eV, the LDOS at the carbon neigh- 
bor of the impurity is broadened and exhibits a peak at 
—8 eV — similar to all covalently bond impurities in their 
minimum energy configuration. The H-impurity causes 
a donor level and is at the same time covalently bond, 
as the supercell band structure with contributions from 
the H impurity marked as fat bands further illustrates. 
There are contributions from the H s-orbital over the 
energy range from -lOeV up to +3eV, indicating strong 
hybridization of the impurity orbital with the graphene 
bands. This is very different from F in the saddle point 
configuration with its valence orbitals contributing sig- 
nificantly within an energy interval which is an order of 
magnitude smaller. In the transition state F and OH are 
ionically bond to graphene. 

The high barrier for H suggests that the formation of 
a strong covalent bond in the transition state is highly 
unfavorable. The origin of this effect can be under- 
stood from the model Hamiltonian, Eqn. [TJ With the 
impurity on top of the bridge site, sublattice symme- 
try is preserved: A = B in Eqn. |T]). The symmet- 
ric combination of the two C-p z orbitals adjacent to 
the impurity <fi + = 775(0, 1, 1) will couple to the impu- 
rity 4>i mp = (1,0,0). The antisymmetric combination 
0_ = -^=(0,1,-1) is decoupled and forms the analog 
of the midgap state occuring for the impurity on top of 
a carbon atom: In the latter case, with the impurity's 
bonding partner in sublattice A, one obtains B = and 
finite A in Eqn. (P). Thus, 4> = (0,0,1) is decoupled 
from the impurity in the stable configuration. (f>o is non- 
bonding and is therefore at the energy of the Dirac point. 
4>- = 773(0, 1' — 1)' however, is an antibonding combina- 
tion of neighboring C-p 2 orbitals. The ab-initio calcula- 
tions show that the resonances derived from this state are 
more than leV above the Dirac point, unoccupied and 
not available for screening the additional positive charge 
brought by the H-impurity. Thus, the creation of a lo- 
cal charge is enforced by graphene's electronic structure 
for the impurity in a b-site saddle point configuration. 
As a consequence, a strong tendency to ionic bonding 
with graphene decreases migration barriers, while migra- 
tion preferably neutral covalently bond impurities is sup- 
pressed. 

This tendency explains experimental findings of 
charged impurities moving almost freely on graphene 



The authors thank A. Geim for inspiring discussions. 
Support from SFB 668 (Germany), FOM (The Nether- 
lands) as well as computer time from HLRN are acknowl- 
edged. 



12, 



29( | and experiments sug gest ing considerable migra- 
tion barriers for H adsorbates [18|, l30( . The fact that clus- 
terization of impurities on graphene strongly suppresses 
their contribution to the resistivity [3l| makes covalently 
bond impurities one natural candidate to main source of 
scattering limiting the electron mobility in graphene. It 
is essential that, as demonstrated here, these impurities 
frequently have quasilocal peaks nearby the neutrality 
point — not accidentally but enforced by symmetry. 
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